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1. The Generalized Eigenvalue Problem 

1.1 History 

At a conference in 1981, K. Wilson suggested to use a variational technique to compute energy 
levels in lattice gauge theory [1]. The idea was picked up and applied to the glueball spectrum [2,3] 
and to the static quark potential(s) [4]. With a certain choice of the variational basis {0,- , / = 1 . . .N} 
and maximizing (0|e^('^''''^|0)/(0|0) with |0) = L,o;,|0,), the variational technique yields the 
generalized eigenvalue problem (GEVP). It is applicable beyond the computation of the ground- 
state energy and has been widely used, but rarely in the form where it can be shown that corrections 
to the true energy levels decrease exponentially for large time [5]. 

Apart from [5], statements about corrections due to higher energy levels seem to be absent in 
the literature. We here add such statements and suggest a somewhat different use of the GEVP, 
which we will show to be more efficient under certain conditions. We also treat the case of an 
effective theory and show numerical results for heavy-quark effective theory (HQET). 

1.2 Basic idea 

We start from a matrix of correlation functions on an infinite-time lattice 

oo 

Cij{t) = {OmOj{t)) = £ e-^"V„;V^„y , i,i = \,...,N (1.1) 

n=l 

Wni = {¥n)i = {n\di\0) En < En+l . 

For simplicity we assume real Yni- States \n) with {m\n) = dmn are eigenstates of the transfer 
matrix and all energies have the vacuum energy subtracted. Oj{t) are any gauge-invariant fields on 
a timeslice t that correspond to Hilbert-space operators Oj whose quantum numbers are then also 
carried by the states \n). Besides the energy levels £"„ one may want to determine a matrix element 

pon = {0\P\n) (1.2) 

of an operator P that may or may not be in the set of operators { } . Starting from the GEVP, 

C{t)Vnit,tQ)=X,^{t,tQ)C{to)Vn{t,to), n = l,...,N t>to, (1.3) 

Luscher and Wolff showed that [5] 

En = limEfit,to), Efit,to) = - log (1.4) 

For a while we now assume that only N states contribute, 

N 



Qjit) = Cfj^ (0 = £ e-^"V„;VA„,- . (1.5) 

We introduce the dual (time-independent) vectors m„, defined by (m,,, V^m) = ^mn ) m,n < N , with 
{un,Ym) = L£i ("«)/ Vmi- Inserting into eq. (U) gives 

d'"(0«„=e-^"V„, d'\t)un=X,['\t,to)d''\to)un. (1.6) 
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So the GEVP is solved by 

?i!l^\t,to)=e-^"^'-'°\ Vn{t,to)ocun (1.7) 

and there is an orthogonality for all t of the form 

(m,„,C(")(Om„) = 5,nnPnit) , Pn{t) = e"^"' . (1.8) 

These equations mean that the operators 2„ = Y!i^i{un)iOi = {0,Un) create the eigenstates \n) = 
Qn\0) of the Hamilton operator: H\n) = En \n) . Consequently we have pon = (0|^|?i) = (0[P2„|0), 
which, preparing for a generalization, we may rewrite as 

V/PMnrm\r ^ L%APit)Ojm)vnit,to)jXn{to + t/2,to) 
j~i (v„(f,fo),C(Ov„(f,fo)) ^ K{tQ + t,tQ) 



while for all t^to we have E^ {t,to) = E, 



n ■ 



Let us now come back to the general case eq. ( |1.1[ ). The idea is to solve the GEVP, eq. (1.3), 



"at large time" where the contribution of states n > N is small and obtain matrix elements and 
energy levels from 



Ef = -log ,y_;'^"|, = E„ + en{t,to) (1.10) 
a An[t + a,to) 

eff i:%i{P{h)Ojm{vn{t,to))jX,{to + t2/2,to) ^ , , n^U 

POn = ; -— -17^ . , , , ^ = POn + nn{t,to) at = = ? • (LH) 

{v„{t,to),C{t2)v„{t,to)r^ Uto + ti,to) 

The restriction to ti = t2 = t is for simplicity. The corrections £„, 7r„ will disappear at large times. 
Note that in the literature the energy levels are often not extracted in this way. Rather, the standard 
effective masses of correlators made from Q„ = {O,v„{t,to)) are used, and the question of the size 
of the corrections is left open. However, the form in eq. ( |1.10| ) has a theoretical advantage as it was 
shown in [5] that (at fixed to) 

e„{t,to) =0{e-'^"'), AE„ = mm\Em-En\. (1.12) 

This is non-trivial as it allows to obtain the excited levels with corrections that vanish in the limit 
of large t, keeping to fixed. However, it appears from this formula that the corrections can be 
very large when there is an energy level close to the desired one. This is the case in interesting 
phenomena such as string breaking [6,7], where in numerical applications the corrections appeared 
to be very small despite the formula above^ . Also in static-light systems the gaps are typically only 
around AE,i ^ 400 MeV, and in full QCD with light quarks a small gap AE^ ~ 2mjc appears in some 
channels. 

Our contribution to the issue is a more complete discussion of the correction e„ to as well 
as a discussion of the corrections n„ to the matrix elements. It turns out that a very useful case is 
to consider the situation 

t<2to, (1.13) 



In fact a different formula was claimed in [6]. 
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e.g. with t — to = const, or 2>t /to = const., and then take to (in practice moderately) large. Then 
it is not difficult to show that 

en{t,to) = O(e-^"+'"0 , A£„,„ =E„-En, (1.14) 
7r„(?,?o) = 0(e"'^"+'"'o), at fixed f-fo (1-15) 
ni{t,to) = o(e-^"+'-''''e-^^2.i(?-^o)) + Q(g-A£„+,.i?) ^ ^^^^^^ 

The large gaps A^T^+i^n can solve the problem of close-by levels for example in the string-breaking 
situation, but also speed up the general convergence very much. For example in static-light systems 
AEd i 2GeV means that roughly a factor of 5 in time separation is gained. We now turn to an 
outline of the proof of these statements. 



2. Perturbation theory 



We start from the solutions above for C = C^^^ and treat the higher states as perturbations. This 
perturbative evaluation was already set up by F. Niedermayer and P. Weisz a while ago [8] but never 
published. We noted the advantage of f < 2to, the form of the corrections to the effective matrix 
elements defined above and could show that these relations hold to all orders in the expansion. 

We want to obtain A„ and v„ in a perturbation theory in e, where 

Avn = XnBvn, A = a'") + eA^'^ , B = B^^^ + eB^^K (2.1) 

We will set 

a(°)=c(°)(0, £AW=cW(0, (2.2) 

=c(°)(fo), eB^'^=d'\to) (2.3) 

in the end. The solutions of the lowest-order equation A^^^v^^ = Xji^^'B^^^ v^"' satisfy an orthogo- 
nality relation (vi^^B^^'vi?) = p„ 5,™ as in eq. (L8) above. Writing 

A„ = xli^^ + eA^^^ + £^Ai^^ . . . , v„ 

we get for the first two orders 



(2.4) 



A(0)v(>)+A(%i°' = Ai 



(0) 



(0) 



+ A„ 



(1) 



b(o),(i)+5(i),(o) 



(2.5) 

+ A,p)B(°)vi°). (2.6) 



With the orthogonality of the lowest-order vectors, v„ , one obtains just like in ordinary QM 

(2.7) 



perturbation theory the solutions for eigenvalues and eigenvectors 



,(1) 



(2) 



-1 -1 
2^ Pn Pm 



(l)^-l/2„(0) 



(1) _ -1/2 
OJfim — Pm 



Vm jA^Vfi 
1 (0) _ . (0) 



I'm ) ''n 

3 (0) _ 3 (0) 
/In yim 



,(0) 



,(0) 



Pn ,i^nVn I I V„ 



,(0) r(i)J0) 



(2.8) 



(2.9) 



Also a recursion formula can be given for the higher-order coefficients. 
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2.1 Application to tlie perturbations C*^') 



Now we insert our specific problem eq. (2.2), eq. (2.3). With straightforward algebra and with 
a representation (for m> n) 

oo 

= (Af )-Hl-e-(^"'-^")('-"'))-i = (Af )-i£e-^(^--^")('-'»), (2.10) 

k=Q 

one finds the correction terms listed at the end of the first section. Initially this is so for the first two 
orders, but the mentioned recursions allow to show that the higher orders are even more suppressed. 

2.2 Effective tlieory to first order 

In an effective theory, all correlation functions 

Qj{t) = C^(t) + (oCjJ'^it) + 0(«2) (2.11) 

are computed in an expansion in a small parameter, (O, which we consider to first order only. The 
notation is taken from HQET where ft) oc l/m. 

We start from the GEVP in the full theory, eq. (1.3), and use the form of the correction terms 
of the effective energies (t < 2to) 

Efit,to) = logj^TT^ = En + 0{c-^^-'-'-% (2.12) 
A„[t + a,to) 

see the discussion above. Expanding this equation in co, we have 

Ef'--{t,to) = a-' log ;^5(" + l"!o) = + O(e-^~?.."0 , (2.13) 

£f .Vm( ) ^ ^^vIt^ - = e!/"" + 0(r e-^^'^..0 . (2.14) 

Here 0{te^^') is a summary for terms {bo + bit)e^^' . As expected for first-order perturbation 
theory, only the eigenvectors of the static GEVP 

a'^\t)vT\t,to) = K''\t,to)C'''\to)vT\t,to), (2.15) 
with normalization (v^^'(f,fo) , C^'^'(fo) vf,'^'(?,?o)) = ^mn , are needed in the formula 

^^"(^^o) = (vr(?,?o), [Ci/'"(0-A,r(f,fo)C'/"^(?o)]vf' (f,fo)) (2.16) 

for the first-order corrections in ft). 
Similarly one may expand 

pf,=pfr' + copf,''/- + Oico^) 

Pof = Pr + 0[e-^""'''"e-^''''"°'(A£yrM^o (2.17) 

and an explicit expression for /?of '^^™ is easily given. Again it involves only the solutions of the 
lowest-order (in (o) GEVP, v^'^' and A,f^', together with the first-order correlators C^/™. The large 
energy gap AE^+i i controls the corrections. 
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3. Application to static-light Bg-mesons 



We have carried out a test in quenched HQET, discretizing the static quark by the HYP2 
action and the strange quark by the non-perturbatively 0(a)-improved Wilson action. Space-time 
is 2L X Lr" with periodic boundary conditions, L k, 1.5fm and we consider two lattice spacings: 
0.1 fm and 0.07 fm (jS = 6.0219 and 6.2885), respectively with k: = 0.133849, 0.1349798. The 
all-to-all strange-quark propagators [9] are constructed from 50 (approximate) low modes and two 
noise fields on each timeslice of 100 configurations. 

The gauge links entering in the interpolating fields are smeared with 3 iterations of (spatial) 
APE smearing [10,1 1]. Then 8 different levels of Gaussian smearing [12] are applied to the strange- 
quark field and we use a simple 70/5 structure in Dirac space for all 8 interpolating fields. The local 
field (no smearing) is included to compute the decay constant. The resulting 8x8 correlation 
function is first truncated to w N y^N one projecting with the N eigenvectors of C(f;) with the 
largest eigenvalues. Here is taken to be roughly 0.2 fm (i.e. = 2a at j8 = 6.0219 and 3a at 
j8 = 6.2885). With N not too large, this avoids numerical instabilities and large statistical errors in 
the GEVP [13]. We present our results for the spectrum and for the decay constant below. 




Figure 1: The estimate aE^ '^''"(f ,?o), n= 1,2, as a function of t, forN = 2,3,4,5 frona top to bottom at a — 0.07 fm. 
The curves are Ef^ + Of^vc '^^'V+i 1 ^ (see comment about AE^-^i i in the text). The coefficients OJ^v are fitted for each N. 



Figure |l] shows the effective energies eq. (1.10) for the lowest two levels at a = 0.07 fm. Sta- 
tistical errors for the ground-state effective energy are below a level of about 3 MeV for time sep- 
arations t < I fm. Unexpectedly, these errors are roughly independent of Iq and of N < 5. The 
functional form of the systematic corrections eq. ( 1.14 ) works very well down to surprisingly small 
t and the independence of to is confirmed by the data. Since the corrections are well understood to 
be below the MeV-level for t > 0.6fm ,N > 4, we may quote for example Ef^'^^ with a total error 
of about 1 MeV. We emphasize that what counts is of course the time separation in physical units. 
The data at the coarser lattice spacing are very similar. 

For this analysis, the energy gaps on the coarser lattice, aAE^+i i ^ 0.46,0.65,0.83, respec- 
tively for N = 2,3,4, have been taken from plateaux of aEf^'^^'^\t,to) for N = 6. They have 
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Figure 2: Bare effective static decay constant as a function of to for different values of t — to at a = 0.07 fm. The curves 
areF + a^e^'^'^+"'° (see comment about AE^^ii inthetext). 



Figure § shows the effective decay constant, eq. ( |1.14 ), at the smaller lattice spacing. The 



leading corrections again dominate at small time already. For N = 5 there is a rather early plateau 
around tQ = 0.4fm, where both excited-state corrections are well below the % level and the sta- 
tistical errors are around 0.7 %. The same statements hold for a = O.lOfm. Note that we fit the 
corrections separately for each t — to and N a.s a. function of tQ. The decay of the fit parameters aM 



as a function oft — to is of the expected form eq. (1.16) 



4. Conclusions 

From a detailed analysis of the corrections to the eigen-values and vectors of the GEVP, it 
becomes clear that tQ should not be made too small. In particular if to > t /2, the simple forms 



eq. ( |1.14| ), eq. ( |1.15| ) can be shown. These corrections decay exponentially with the large gaps 
En^\ — En. For first-order corrections in an effective theory a similar suppression holds, with the 
energy differences of the lowest-order theory. 

As pointed out to us at the conference, the authors of [14] studied the GEVP for a toy model 
with ten states and noted that it is relevant to have tQ "large enough". Fig. 17 of [14] indeed illus- 
trates that the effective energies become independent of tQ when (roughly) tQ>t/2is respected. 
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